HAMILTONIAN DESCRIPTION AND TRAVELING WAVES OF THE 

SPATIAL DYSTHE EQUATIONS 



FRANCESCO FEDELE* AND DENYS DUTYKH 



Abstract. The spatial version of the fourth-order Dysthe equations describe the evolu- 
tion of weakly nonlinear narrowband wave trains in deep waters. For unidirectional waves, 
the hidden Hamiltonian structure and new invariants are unveiled by means of a gauge 
transformation to a new canonical form of the evolution equations. A highly accurate 
Fourier-type spectral scheme is developed to solve for the equations and validate the new 
conservation laws, which are satisfied up to machine precision. Further, traveling waves 
are numerically investigated using the Petviashvili method. It is found that their collision 
appears inelastic, suggesting the non-integrability of the Dysthe equations. 
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1. Introduction 

Mathematical models used in physics and mechanics do not always possess a canonical 
Hamiltonian structure. Typically, the dynamics is governed by partial differential equa- 
tions expressed in terms of physically-based variables, which are not usually canonical. A 
transformation to new variables is needed in order to unveil the desired structure explicitly 
(see, for example, [24]). This is the case for the equations of motion for an ideal fluid: in the 
Eulerian description, they cannot be recast in a canonical form, whereas in a Lagrangian 
frame the Hamiltonian structure is revealed by Clebsch potentials (see, for example, [24], 



Key words and phrases. Dysthe equation; deep water; Hamiltonian; ground state; traveling wave; soli- 
tary wave. 

* Corresponding author. 

1 



2 



F. FEDELE AND D. DUTYKH 



[21]). Moreover, multiple-scale perturbations of differential equations expressed in terms 
of non-canonical variables typically lead to approximate equations that do not maintain 
the fundamental conserved quantities, as the hydrostatic primitive equations on the sphere 
[19], where energy and angular momentum conservation are lost under the hydrostatic 
approximation. 

Clearly, if canonical variables can be identified, then the associated Hamiltonian struc- 
ture provides a natural framework for making consistent approximations that preserve the 
fundamental dynamical properties of the original system, notably its conservation laws 
[23, 25]. For example, consider the equations that describe the irrotational flow of an ideal 
incompressible fluid of infinite depth with a free surface. Their Hamiltonian description 
was discovered by [35] in terms of the free-surface elevation r)(x, t) and the velocity poten- 
tial <p(x, t) = (f>(x, z = r)(x, t),t) evaluated at the free surface of the fluid. Variables r)(x, t) 
and <p(x,t) are conjugated canonical variables with respect to the Hamiltonian % given 
by the total wave energy. By means of a third order expansion of % in the wave steep- 
ness, [36] derived an integro-differential equation in terms of canonical conjugate Fourier 
amplitudes, which is Hamiltonian and has no restrictions on the spectral bandwidth. To 
derive the Zakharov equation, fast non-resonant interactions are eliminated via a canonical 
transformation that preserves the Hamiltonian structure [16, 36]. Alternatively, envelope 
equations for the free-surface and wave potential can be obtained via multiple-scale per- 
tubation techniques that remove such terms [35, 28]. However, the Hamiltonian structure 
is lost since the transformation from the Zakharov action variables to envelope ampli- 
tudes is not canonical. Nevertheless, the Hamiltonian and non-Hamiltonian versions of the 
Zakharov's equation are completely equivalent from a practical view-point. 

The modified Nonlinear Schrodinger (NLS) equations derived by Dysthe [7] are also non- 
Hamiltonian. Using the method of multiple scales, he extended the deep-water cubic NLS 
equation for the time evolution of the unidirectional narrowband envelope A of the velocity 
potential with carrier wave exp(ikox — iuot) and that of the potential <fi of the wave-induced 
mean flow, to fourth order in steepness and bandwidth. Introducing dimensionless units, 
t' = uot, x' = koX, A' = k A, and dropping the primes, (ft can be easily found by means of 
the Fourier-transform and a single equation for A can be derived as well: 

A t + \a x + l -iA xx - ^A xxx + % - \Af A^-\A\ A x + (3A 2 A* X + UaH(\A\ 2 ) x = 0, (1.1) 

where /3 = 1/4, the subscripts A t = d t A and A x = d x A denote partial derivatives with 
respect to x and t respectively, H(/) is the Hilbert Transform of a function f(x), and A* 
denotes complex conjugation (see also [14]). The form of the equation for the envelope B 
of the free surface is similar to (1.1), but the term (3A 2 A* becomes —(3B 2 B* [12, 13]. As 
noted by [12] such terms do not play any role in the linear stability of unifrom wave trains 
within the order of truncation of the perturbation analysis leading to the Dysthe equation. 
However, in three dimensions the /3-term yields a better confinement of the instability 
region for the Stokes wave and the associated Dysthe equations are less susceptible to 
numerical instability and energy leakage [31]. Moreover, [13] noted that the two equations 
for A and B are so similar 'because of some deeper structure in the problem which is as 
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yet unresolved 1 . This may well be a Hamiltonian structure that the temporal Dysthe 
equations lack because (1.1) is expressed in terms of envelope variables, which are not 
canonical. Recently, the associated canonical form (which do not contain the /3 term) has 
been derived by [11] starting from the Hamiltonian Zakharov equation with the Krasitskii 
kernel [16]. We also point out that [37], starting from a conformal-mapping formulation 
of the Euler equations derived another version of the temporal Dysthe equation, which is 
similar to (1.1) but also non-Hamiltonian. 

On the other hand, to model wave propagation in wave basins a change to a coordinate 
system moving at the group velocity can be used by introducing the dimensionless variables 

A = ev, r = e(2x — t), £ = e 2 x, 

with e = k a being the wave steepness of the carrier wave and a the associated amplitude 
[18]. As such, the temporal Dysthe equation (1.1) transforms, up to the fourth order in e, 
to 

i>£ + iv TT + i \v\ 2 v + 8e \v\ 2 v T + 2ewH(\v\ 2 ) T = 0. (1.2) 

hereafter referred to as the spatial Dysthe of the envelope v of the wave potential. On the 
other hand, the associated envelope B = eu of the free surface satisfies (see also [31]) 

+ iu TT + i \u\ 2 u + 8e \u\ 2 u T + 2eu 2 u* T + 2eiuH(\uf) T = 0, (1.3) 

Both (1.2) and (1.3) can also be derived directly from the Zakharov equation [15]. 

In this paper, we will unveil the hidden canonical structure of (1.2) and that of (1.3). 
Following our previous letter [8] , we first review the properties of Hamiltonian systems and 
then introduce a gauge transformation that yields a canonical form for the spatial Dysthe 
(1.3), and new invariants for it. As a corollary, we will also show that equation (1.2) for the 
wave potential envelope v is already Hamiltonian, as expected since (1.2) is a special case 
of the more general result obtained for two horizontal dimensions by [4]. Further, a highly 
accurate pseudo spectral scheme is exploited to validate the new conservation laws up to 
machine precision. Then, ground states and traveling waves are numerically investigated 
using the Petviashvili method ([22]; see also [17], and [34]). Their collision is numerically 
investigated in order to provide new insights on the integrability of the spatial Dysthe 
equations. 



2. Hamiltonian description 

For the spatial evolution of canonical systems, the Hamilton's equations take the form 
(see, for example, [3, 21]) 

where u l (£, r) is the i th variable of the dynamical vector u, = d^u 1 denotes space 
derivative, "H(u, £) is the Hamiltonian, and J l J is a skew-symmetric (symplectic) matrix. 
The Einstein's notation for repeated index summation is adopted, and 6 denotes variational 
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different at ion. If % does not depend explicitly on £, the flow generated in the phase-space 
of u conserves %. This follows from (2.1) 

— = d H + —d u j = ™Lfi 5 -H = o 

d£ ^ 5uj ^ 5uj c 5uj ' 

since d(H. = and J l J = —J ] c l . If one considers a general transformation 

w i = u/(u), (2.2) 
the Hamilton equations (2.1) in the new variables become 

<W = J im (w)-^-, (2.3) 

where 

dw l dw 3 

v ' du n c du m y ' 

is the cosymplectic form (see, for example, [21]). The transformation (2.2) is canonical if 

jv = ^Ljnm^L (2.5) 
du n c du m v ' 

then the Hamilton's equations (2.1) are invariant, otherwise (2.3) is in a noncanonical 
form and J'- 7 may be singular. However, if it satisfies certain algebraic properties, i.e. 
skew symmetry and the Jacobi identity, then a generalization of the Darboux theorem 
garanties the existence of a transformation from J u to J l J (see, for example, [21]) for 
a more comprehensive discussion on canonical and noncanonical forms). Obviously, by 
switching £ with r the above description also holds for the time evolution of Hamiltonian 
systems. 

As an example, consider the finite dimensional system of an harmonic oscillator in the 
classical canonical variables q(t) (coordinate) and p(t) (momentum). This admits the 
canonical form 

. &h m 

Q = ^-=P, P = — 5- = -?» ( 2 - 6 ) 
op oq 

where p denotes time derivative, and Ti = (q 2 + p 2 )/2. The flow in the phase-space is 
'incompressible' since the divergence vanishes: 

dq dp 

— + — = 0. 

dq dp 

The transformation z = q + ip is canonical and (2.6) transforms to 

, = -z— = - % z, z = z— = iZ , (2.7) 

where Ji = \z\ 2 /2, and z* is the complex conjugate of z. It is straightforward to prove 
that the gauge transformation 

z = we iaH \ 
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with a as a free parameter, is also canonical, and (2.7) remains unchanged in the new 
variables w and w*. On the other hand, if one considers the coordinate change 

Q= 7== . P = V, (2.8) 



then (2.6) transforms to the noncanonical form 

P = -Q (Vl + aP* - 2aP 2 ) , Q= J—- . (2.9) 
^ ' v 1 + olP 

This flow does not preserve volume as (2.6) does, nonetheless equations (2.9) are those 
of a disguised harmonic oscillator obtained via the noncanonical change of variables (2.8) 
(see also [21]). The Dysthe equation (1.3) shares the same roots as (2.9). They both come 
from a noncanonical transformation of a Hamiltonian system. In the following, canonical 
variables for the spatial Dysthe are unveiled. 

3. Canonical form of the spatial Dysthe equations 

Hereafter, we will consider the generic nonlinear equation 

= —iau TT — ih |w| 2 u — ce \u\ 2 u T — eeu 2 u* T — fieuH(\u\ 2 ) T , (3.1) 

with (a, h, c, e, /) as a quintuplet of arbitrary real coefficients. In particular, the spatial 
Dysthe for the wave and potential envelopes follow from (3.1) with parameters (1, 1, 8, 2, 2) 
and (1, 1, 8, 0, 2) respectively. The wave action 



A= I \u\ 2 dr (3.2) 



is conserved by (3.1) since 
where 



d{ \uf + d T R = 0, (3.3) 

K = ia[u T u — u T u) H — e \u\ . 

Up to date, no other conservation laws are known for u. Drawing from [6] (see also [33]), 
the invariance of A suggests the following variable change via the gauge transformation 

w = Q{u) = u exp(ikip), (3.4) 

where k is a free parameter, and the 'stream function' ip is defined as 

Note that \u\ 2 = \w\ 2 and the wave action A is preserved in the transformation. Further, 
from (3.3) the spatial evolution of ^ is governed by 

d^ T ip = —d T R, 

so that dgij) = ?/>g = —R. The spatial evolution equation for w follows from (3.4) as 

W£ = exp(ihip) + ikwip£ = exp(ikip) — ikwR. (3.5) 
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From (3.1), can be given in terms of time derivatives of u, which by means of (3.4) can 
be expressed in terms of w and ip as 

u T = exp(-ikip) (w T — ikwip T ) , 

u TT = exp(— 2/V?/>) (w TT — ikwip TT — 2ikw T ip T — k 2 wipl) . 
Explicitly, the evolution equation for w is given by 

. . 2 7 2^ 2dk . .4 

W£ = —iaw TT — ih\w\ w + ike \w\ w 

— (c + 2ak)e \w\ w T — (e + 2ak)ew 2 w* — ifewH.(\w\ 2 ) T . (3.6) 
If the free parameter k* is chosen as 

equation (3.6) simplifies to 

= —iaw TT — ih \w\ 2 w — 2— — ^ — — £ 2 \w\ 4 w — (c — e)e \w\ 2 w T — ifewH(\w\ 2 ) T , (3.7) 



which admits the Hamiltonian structure 



W H \ _ ■ ( 1 



wt V —1 



Jiw 

? / V /V <5ui* 

where the Hamiltonian is given by 

2 h . .4 ce — 2e 2 2 , , 
- - \w\ — — e \w\ 



(3.1 



i — - — e |w| 2 (w*w — w T w*) — — e \w\ 2 H(|w| 2 ) T )dr. (3.9) 



4 1 1 v r ' J 2 
Here, % w is also an invariant together with the momentum 



M w = J i(w*w - w T w*)dr. (3.10) 

Note that in the Appendix of [31] it was shown that Eq. 1.1 admits the same invariant 
(3.10) if (3 = 0. Since H w is an invariant for w, from (3.4) and (3.10) 

S{u) = H w (g-\w))= f(a\u T \ 2 -^\u\ 4 + ^±^e 2 \u\ 6 (3.11) 

J V 2 6a 

■° + e I |2 / * *\ / 1 |2tt/i i2\ \ j 

—2 — - — £ |«| (M T M — M T M J — — £ |M| H(|lt| ) T jdT, 



M{u) = M w {Q-\w)) 



g 

i(u*u — u T u*) e \u\ 

a 



are both invariants of the original spatial Dysthe equation (3.1). Unfortunately, £ is not 
the associated Hamiltonian. In section 5 we present a highly accurate Fourier-type pseudo- 
spectral method that has been used to solve for the envelope dynamics and validate the 
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invariance of the Hamiltonian (3.9) of w and the new invariants (3.11) of u. In our numerical 
investigations, it is found that they are all conserved up to machine precision. 

Note the appearance of terms of 0(e 2 ) in the invariants £ and % w . They both vanish if e 
is null; as a result u = w and £{u) = H w becomes the Hamiltonian for u. As a consequence, 
the spatial Dysthe (1.2) for the wave envelope B is non Hamiltonian since e = 2, whereas 
that for the wave potential A is (cf. Eq. 1.3)). 

To the leading order, the free-surface rj is related to the wave envelope A = eu as 

k rj ~ eue i{kox -" ot \ 
where from (3.4), u can be expressed in terms of the canonical w as 

k T] ~ £we i{k^- Uo t-iek*\i>\) _ (3 12 ) 

Thus, via the gauge transformation (3.4) the wavenumber shift due to self-resonant interac- 
tions naturally arises, and the envelopes \u\ and \w\ are identical. This may suggest (3.12) 
as a canonical ansatz for multiple scale perturbations of the Euler equations. Further, w 
satisfies the canonical equation (3.7), viz. 

= —iw TT — i \ w\ 2 w — 2ie 2 \w\ 4: w — 6e\w\ 2 w T — 2iewH(\w\ 2 ) T , (3.13) 

with the Hamiltonian given by 

n I I ( I |2 1 I |4 • 3 I |2/ * *\ 

H w — I \\w T \ — -\w\ — i—e\w\ [w T w — w T w ) 

-eM 2 H(M 2 ) T - ^ 2 |ufW (3.14) 

which is invariant together with the action (3.2) and momentum (3.10). Note that the 
corresponding 0(e 2 ) terms in (3.13) and (3.14) can be legitimately neglected and the 
Hamiltonian properties are preserved [23, 25]. 

The noncanonical envelope u satisfies, besides the action (3.2), the two new invariants 

£{u) = J ^\u T \ 2 — — \u\ A — i^e \u\ 2 (u*u — u T u*) + ~^ £2 \ u f ~ e l u | 2 -^(l n | 2 )r^ dr, 

M{u) = J [i(u*u — u T u*) — 2e |m| 4 ] dr, 

Practically, terms of 0(e 2 ) can be neglected and the invariants are approximately satisfied. 
The spatial Dysthe (1.3) of the wave potential A = ev is already Hamiltonian with 

H„ = y (^\v T \ 2 -^\v\ 4 -2ie\v\ 2 (v*v -v T v*) - e\v\ 2 H(\v\ 2 ) T ^J dr, 

since the term e = 0. Besides / H V , it also conserves the action (3.2) and momentum (3.10) 
written in terms of v. 



8 



F. FEDELE AND D. DUTYKH 



4. Ground states and Traveling waves 



Insights into the underlying dynamics of the Dysthe equations are to be gained if we 
construct some special families of solutions in the form of ground states and traveling 
waves, often just called solitons or solitary waves. Hereafter, we do so for the equation 
(3.1) for the envelope u, and the associated Hamiltonian form (3.7) in w can be treated in 
a similar way. However, owing to the gauge transformation (3.4) the envelope \w\ = \u\. 

First, consider the special case of (3.1) with parameters (1, h, 0, 0, 0). This is the classical 
NLS equation, which is well known to admit a family of localized traveling waves of the 
form ([38]; see also [1, 29]) 



where A is a free parameter. This solitary wave travels along the space £ with speed 1/s 
while pulsating in time. Due to the integrability of the NLS equation the interaction of 
solitary waves is elastic and so after the collision their shape is unchanged, but a phase 
shift occurs (see, for example, Figure 5). At infinite speed, i.e. s = 0, (4.1) yields the 
so-called ground state solutions (see, for example, [34]). 

In the following, we wish to find traveling waves of the Dysthe equation (3.1) in u. To 
do so, we consider a generalization of (4.1) via the ansatz 



where /x and s are generic parameters and the function v(-) is in general complex. Substi- 
tuting (4.2) into (3.1) and separating the real and imaginary parts v = $ + yield 

/x$ - a$ TT + s^/ T = /i$($ 2 + ^ 2 ) + c^ T ($ 2 + ^ 2 ) 

+e(2$$ r ^ - $ 2 ^ r + tf 2 tf T ) + /$H(($ 2 + ^ 2 ) r ), (4.3) 
/i*-a* Tr -s$ T = /i*($ 2 + ^ 2 ) -c$ T ($ 2 + ^ 2 ) 

-e($ 2 $ r - $ T ^ 2 + 2$^^ T ) + /^H(($ 2 + ^ 2 ) T ). (4.4) 



These equations can be solved analytically if the c, e and / terms can be assumed as small 
perturbations of the NLS equation. In this particular case direct soliton perturbation 
theory can be applied (see, for example, [2] and [34]). We prefer to deal with the general 
case, and thus numerically solve for (4.3) and (4.4) using the Petviashvili method ([22], see 
also [17, 34]). This numerical approach has been successfully applied by [37] to compute 
ground states of their version of the temporal Dysthe equation. 

To shorten the notation we rewrite the differential equations (4.3) and (4.4) in the 
operator form separating linear and nonlinear terms: 




(4.1) 



u(C, t) = v(t - s£)e 



(4.2) 




SPATIAL HAMILTONIAN DYSTHE EQUATION 



9 



where J\f(Q, \P) denotes the right-hand side of equations (4.3), (4.4) and the matrix C is 
defined as 



To initialize the iterative process, one can use the analytical solution to the NLS (4.1) 
for the real part $o and taking the imaginary part ^0 = 0. This iterative procedure 
is continued until the L OT norm between two successive iteration is of the order of the 
machine precision. We would like to underline the fact that this method can be very 
efficiently implemented using the Fast Fourier Transform (FFT) (see, for example, [9]). 

As an application, consider the non-Hamiltonian Dysthe (1.2), particular case of (3.1) 
with parameters (1, 1, 8, 2, 2). Figure 1 shows the action A of the ground state (s = 0) for 
u (and so w due to the gauge transformation (3.4)) as function of fi, computed for different 
values of the steepness e, and Figure 2 reports the associated envelopes \u\ = \w\ for /i = 3. 
As one can see, as e increases they tend to reduce in size in agreement with the asymptotic 
analysis carried out by [2] in the limit of e — > 0. 

Furthermore, Figure 3 illustrates a typical basin of attraction of the Petviashvili scheme 
in the phase space (s,fi) for e = 0.15. Each gray dot corresponds to a well converged 
solution up to machine precision, whereas white spots are associated to either divergent or 
converged-to-zero solutions. In particular, we noted that the numerical scheme converged 
to localized traveling waves below the black boundary curve T shown in Figure 3. On the 
other hand, the convergence to simple periodic waves occurred for points above T. An 
analytical form of such curve can be derived from the nonlinear dispersion relation A - u 
of periodic waves of the form 




where C is the symbol of the operator C and k is the Fourier transform parameter. Then 
the Petviashvili iteration takes the following form: 




operator £ 1 in the Fourier space is given by the following matrix: 




u(£,t) = Ae^-^e-^. 
Indeed, using this ansatz in (3.1) yields 

au 2 + (s + (e - c)A 2 )u + /i - hA 2 = 0. 



(4.5) 
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Figure 1. Action .A of the ground state of the spatial Dysthe equation (1.2) 
for u as function of /i, for different values of the steepness e. 

The boundary T of the periodic waves region in the phase space (s,/i) follows from the 
condition that the discriminant of the quadratic equation in u is null, viz. 



The analytical separatix (4.6) is practically identical to the boundary T that we have 
identified numerically, confirming the validity and correctness of our numerical solver. 
Thus, localized traveling waves bifurcate from T. This is clearly illustrated in Figure 4, 
which reports the change in shape of the envelope |m| as s varies while keeping \i constant, 
say fi = 0.27 as a particular case. Clearly, as the point (s,fi = 0.27) reaches the boundary 
T from below (see Figure 3), the soliton envelope tends to flatten to that of a periodic 
wave. Similar results also hold for the Hamiltonian w due to the gauge transformation 



We further investigate the dynamics of the Dysthe equations by means of a highly ac- 
curate Fourier-type pseudo-spectral method, which is needed in order validate the newly 
derived conservation laws and investigate nonlinear phenomena such as elastic/inelastic 
traveling waves interactions. Spectral methods have been already proven to be sufficiently 
accurate and robust for studying the dynamics of gravity waves [20, 5, 34]. As in the previ- 
ous section, we first illustrate the numerical method for equation (3.1), but the associated 




(4.6) 



(3.4). 



5. Traveling wave interactions 
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Figure 2. Envelopes \u\ of ground states of the spatial Dysthe equation 
(1.2) for different values of the steepness s (fi = 3). The associated canonical 
envelope \w\ = \u\ because of the gauge transformation (3.4). 

Hamiltonian form (3.7) can be treated in a similar way. Then we present some results of 
the interactions among traveling waves. 

5.1. Numerical method description. Equation (3.1) can be written in the following 
operator form: 

u^ + C-u =N(u), (5.1) 

where the operators L and H are defined as 

£ = iad TT , Af(u) = — ih\u\ 2 u — ce\u\ 2 u r — eeu 2 u* — ifeuH(\u\ T ). 

Equation (5.1) will be solved numerically by applying the Fourier transform in the time 
variable r. The transformed variable will be denoted by u = J-{u}. We recall that the 
symbol of the last term H(<9 r ) is equal to \k\, and that of L is —iak 2 , k being the Fourier 
transform parameter. The nonlinear terms are computed in the physical space, while time 
derivatives and the Hilbert transform are computed in the Fourier space. For example, the 
term |n| 2 n T is discretized as: 

F{\u\ 2 u T } = ^{l^" 1 ^)! 2 • F^iiku}}. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 



V- 

Figure 3. Numerical basin of attraction of the Petviashvili scheme in the 
phase space (s,/x) for e = 0.15. Each gray dot corresponds to a well con- 
verged solution up to machine precision, whereas white spots are associated 
to divergent or converged-to-zero solutions. Solitary waves (localized trav- 
eling waves) occur below the black boundary curve T, which separates the 
region of periodic waves. 

The other nonlinear terms are treated in a similar way. We note that the usual 4-half rule 
is applied for anti-aliasing [30, 5, 10]. 

In order to improve the stability of the space discretization procedure, we can integrate 
exactly the linear terms. This is achieved by making a change of variables [20, 10]: 

Finally, the resulting system of ODEs is discretized in space by the Verner's embedded 
adaptive 9(8) Runge-Kutta scheme [32]. The step size is chosen adaptively using the so- 
called H211b digital filter [26, 27] to meet the prescribed error tolerance, set as of the order 
of machine precision. 

5.2. Numerical results. As an application, consider the interaction between two solitary 
waves of the NLS equation traveling in opposing directions with the same speed s = 2, for 
fi = 2 (e = 0). The shape of the two solitons is identical because of the NLS reflection 
symmetry, i.e. r) = (£, — r). The left plot of Figure 5 shows that the two solitons 
emerge out of the collision with the same shape, but a phase shift. The interaction is 
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Figure 4. Localized traveling waves bifurcating from the boundary V of 
Figure 3. Change in shape of the envelope \u\ as s varies while keeping 
/i = 0.27. In particular, with reference to Figure 3, the point (s,/i = 0.27) 
reaches the boundary V from below. 



elastic as it should be since the NLS equation is integrable. This is clearly seen from 
the right plot of the same figure, which reports the initial and final shapes of one of the 
two solitary waves. We point out that the Hamiltonian and the other two NLS invariants 
are conserved up to machine precision. For the same parameters, Figure 6 reports the 
interaction of the associated solitary waves of the Dysthe equation (1.2) for e = 0.15. The 
reflection symmetry is lost and the two solitons have different shape and amplitude. The 
interaction is clearly inelastic, since after the collision radiation is shed and the initial and 
final soliton shapes are different as seen in the right plot of Figure 6. The interaction of 
four solitons has similar inelastic characteristics as shown in Figure 7. This suggests the 
non-integrability of the Dysthe equation (1.2). Similar dynamics is also observed for the 
associated Hamiltonian form (3.7). 
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x 



Figure 5. (Left) elastic collision of two NLS solitary waves traveling at the 

same speed s — 2, for \i — 2, and (right) initial ( ) and final ( — ) shapes 

of one of the two solitons. 



6. Conclusions 

A canonical variable for the spatial Dysthe equation for the wave envelope u has been 
identified by means of the gauge transformation (3.4), and the hidden Hamiltonian struc- 
ture is unveiled. Further, the associated equation for the wave potential v is shown to be 
already Hamiltonian. Exploiting the gauge invariance allowed to derive two new invariants 
for the noncanonical u. Clearly (3.4) is also applicable to the temporal Dysthe equation 
(1.1), but it leads to a cumbersome equation in w that contains nonlinear terms with 
derivatives higher than the first order, which cannot be eliminated by a proper choice of 
the free parameter k. The transformed equation may well be Hamiltonian, but we could 
not identify such functional due to the complexity of the equation. In any case, the gauge 
transformation (3.4) does not yield the Hamiltonian form derived by [11] for the temporal 
case. 

Moreover, the existence of localized traveling waves that bifurcate from periodic waves 
has been investigated numerically by means of a highly accurate Petviashvili scheme. In 
particular, ground state solutions are found to be in agreement with the asymptotic analysis 
of [2] . A Fourier-type pseudo-spectral method has been developed to solve for the envelope 
dynamics and test the new conservation laws, which are validated up to machine precision. 
Solitary waves are found to interact inelastically, suggesting the non-integrability of both 
the Hamiltonian and non-Hamiltonian version of the spatial Dysthe equations. 
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Figure 7. Inelastic collision of four Dysthe solitary waves. 
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